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Abstract 

Exact solutions of the equation governing the equilibrium magnetohydrody- 
namic states of an axisymmetric plasma with incompressible flows of arbitrary 
direction [H. Tasso and G.N. Throumoulopoulos, Phys. Pasmas 5, 2378 (1998)] 
are constructed for toroidal current density profiles peaked on the magnetic axis 
in connection with the ansatz S = —ku, where S = d / du [g(d<& / du) 2 ] (k is a 
parameter, u labels the magnetic surfaces; q(u) and are the density and the 
electrostatic potential, respectively). They pertain to either unbounded plasmas 
of astrophysical concern or bounded plasmas of arbitrary aspect ratio. For k — 0, 
a case which includes flows parallel to the magnetic field, the solutions are ex- 
pressed in terms of Kummer functions while for k 7^ in terms of Airy functions. 
On the basis of a tokamak solution with k 7^ describing a plasma surrounded 
by a perfectly conducted boundary of rectangular cross-section it turns out that 
the Shafranov shift is a decreasing function which can vanish for a positive value 
of k. This value is larger the smaller the aspect ratio of the configuration. 
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1. Introduction 



Flow is a common phenomenon in astrophysical plasmas, in particular let us 
mention here the plasma jets that are observed to be spat out from certain galactic 
nuclei and propagate to supergalactic scales, e.g. Ref. 0. Also, there has been 
established in fusion devices that sheared flow can reduce turbulence either in the 
edge region (L-H transition) or in the central region (internal transport barriers) 
thus resulting in a reduction of the outward particle and energy transport, e.g. 
Ref. 0. 

In an attempt to contribute to the understanding of the equilibrium proper- 
ties of flowing plasmas we considered cylindrical axisymmetric jSlEl) an d 
helically symmetric j7j steady states with incompressible flows in the framework 
of ideal magnetohydrodynamics (MHD) by including the convective flow term in 
the momentum equation. Also we studied the equilibrium of a gravitating plasma 
with incompressible flow confined by a point-dipole magnetic field [H] . For an ax- 
isymmetric magnetically confined plasma the equilibrium satisfies [in cylindrical 
coordinates (z, R, <fi) and convenient units] the elliptic differential equation for 
the poloidal magnetic flux-function ip, 

(l-M p ^-i K)<|V^ + i (jiy + rf = 0, (1) 

along with a Bernoulli relation for the pressure, 



PsW - Q 



v 



2 t?2frf>'\2 



(2) 



Here, P s (?p), p(ip) and $(/0) are, respectively, the static-equilibrium pressure, 
density and electrostatic potential which remain constant on magnetic surfaces 
i/j(R, z) = constant; the flux functions F(if)) and X(ip) are related to the poloidal 
flow and the toroidal magnetic field; M p = F'/y/g is the Mach-number of the 
poloidal velocity with respect to the poloidal-magnetic-field Alfven velocity; A* = 
i? 2 V • (V/i? 2 ); the prime denotes differentiation with respect to ip. For vanishing 
flow ((TJ) and (J2J) reduce to the Grad-Schliiter-Shafranov equation and P = P s (ip), 
respectively. Derivation of (0) and (j2J is given in Ref. [S]. It should be clarified 
that to simplify (22) of [5], the flux-function term P s - XF'<&'/(1 - M 2 ) therein 
has been replaced by P s in (1) here. Consequently, relation (19) of [5] for the 
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pressure takes the form (2) here. Also, to emphasise that a poloidal-velocity Mach 
number is involved, M in [5] is denoted by M p here. 
Under the transformation 



(JTJ) reduces to 



u(ij>) = £ [1 - Mp(g)) l/2 dg, M 2 p < 1, 



2du \ 1 - Ml ) + ~d^ + Ydu 



(3) 




(4) 



The flow contributions here are connected with Mp and S 



d/du g(d<&/du] 

Equation (j3J), free of the nonlinear term l/2(M p 2 )'\ V?/>| 2 , can be analytically solved 
by assigning the ^-dependence of the free flux functions P s , X, F, g, and 
Relation (J2J) then determines the pressure. For S = constant, exact solutions 
which extend the well known Solovev one were derived in Ref. jH]. Owing to the 
flow and its shear a variety of new configurations are possible having either one or 
two stagnation points in addition to the usual ones with a single magnetic axis. 

The aim of the present work is to construct exact solutions of (jlj) for S oc u 
and examine their properties. The solutions are of the form (JJJ) below in Section 
2. This form is advantageous in that boundary conditions associated with ei- 
ther bounded laboratory plasmas or unbounded plasmas of astrophysical concern 
can be treated in a unified manner by adjusting appropriately the parameters 
it contains. As a matter of fact certain of the solutions we construct constitute 
extensions of the Hernegger-Maschke solution of the Grad-Schluter-Shafranov 
equation for bounded plasmas [3 HO] and a recent one for unbounded plasmas 
derived by Bogoyavlenskij [TTj . 

The work is organized as follows. Exact solutions of (jlj for S oc u are con- 
structed in Section 2. The flow impact on the Shafranov shift is then examined in 
Section 3 on the basis of a particular solution describing a tokamak configuration 
of arbitrary aspect ratio being contained within a perfectly conducting boundary 
of rectangular cross-section. Section 4 summarizes the Conclusions. 
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2. Exact equilibrium solutions 



Let us make the following particular choice for the free functions: 

X 2 / 

X 2 + cju 2 , P s = P s0 + 2(-l) 7 c 2 ! w 2 , and S = p($') 2 = -QAc 3 c 3 2 u 



1-M 2 



(5) 

under which (J3J) becomes linear. Here, the factor (— l) 7 has been introduced in 
order to make comparison with solutions existing in the literature convenient, 
and X , P s0 , c 1; c 2 , c 3 are constants. Equation (jlj) then takes the form 



A*u + c\u + 2(-1) 7 c 2 xm - 8c 2 c 3 x 2 u = 0, 
where x = 2c 2 R 2 . We pursue separable solutions of the form 

u(R, z) = x n P{x)T{z) exp (—) ■ 



(6) 



(7) 



This is appropriate for considering various equilibrium configurations in connec- 
tion with different boundary conditions. In particular, the term x n which makes u 
to vanish on the axis of symmetry is associated with either compact tori or spher- 
ical tokamak configurations. Plasmas surrounded by a fixed perfectly conducting 
boundary can be considered by setting 7 = 0. Unbounded plasmas are connected 
with 7 = 1, the exponential term guaranteeing smooth behaviour at large dis- 
tances. It is noted that for static equilibria (M p = C3 = 0) the cases n = 7 = 
and (n = 0,7 = 1) were considered in Refs. jHHD], and [TT], respectively. 

With the use of (J7J) equation (jHJ) leads to the following differential equations 
for P(x) and T(z): 

T" + u 2 T = (8) 

and 



+ 



4n(n — 1) — 47nx + 7 2 x 2 



+ 



Ax 4 
where uj 2 = constant and r = (to 2 — c 2 )/8c2. Therefore, 



xP" + (2n - 7 x)P' 
P = 0, 



l) 7 x 2 

- T - C 3 X 



(9) 



T(z) = A sm(ujz) + B cosfu;;?) 



(10) 
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and the problem additionally requires to solving Q. Solutions of (jHJ) associated 
with different boundary conditions will be constructed in the subsequent subsec- 
tions. It is noted here that all the solutions hold for arbitrary poloidal Mach 
numbers, viz. the dependence of M 2 {u) on u remains free. 

2.1 Unbounded plasmas (j = 1) and C3 = 

It is first noted that the case C3 = includes flows parallel to the magnetic 
field for $' = 0, viz. when the electric field vanishes. For C3 = 0. (fl|) becomes 
identical in form with the Grad-Schliiter-Shafranov equation and Q reduces to 

x 2 P" + x{2n - x)P' + [n(n - 1) - (r + n)x\ P = 0. (11) 

Equation (fTTj) can be solved by the following procedure. The substitution P = 
x k W(x), where k is a root of the quadratic equation 

k 2 + (2n- l)k + n(n- 1) = 0, (12) 

leads to the following equation for W(x): 

xW" + [2(Jfe + n) - x}W -(k + r + n)W = 0. (13) 

The solutions of this equation can be expressed in terms of Kummer or confluent 
hypergeometric functions ([12!, P- 503; [13], p. 137). In particular, for the two 
roots of (fT2|) we have: 
2.1. a ki = —n 

The solution of (|T3|) is 

W(x) = x [DiM(t + 1, 2, x) + D 2 U(t + 1, 2, x)] , 

where M and U are the Kummer functions of first and second kind, respectively, 
and Z?i, D 2 are constants. Consequently, the solution of the original equation (JJJ) 
is written in the form 

u(x,z) = x[D 1 M(T + l,2,x) + D 2 U(T+l,2,x)]exp(-x/2) 

[A sin(ujz) + B cos(luz)} . (14) 

For special values of r and n the Kummer functions reduce to simpler classical 
functions or polynomials (see Ref. |T2|, p. 509, table 13.6) ; in particular, for 
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t = —m, where m is a non-negative integer, and n = they reduce to Laguerre 
polynomials. Solutions of this kind for static equilibria (M p = 0) were derived 
in Ref. ^T] and employed to model astrophysical jets and solar prominences. It 
may also be noted here that a continuation of this study and related studies were 
reported in Refs. [HJ[T3] and [15] . 

2.1.1 k 2 = -(n + l) 

The solutions of (fT3J) and (jH), respectively, read 

W(x) = x 3 [D 1 M(t + 2, 4, x) + D 2 U(t + 2, 4, x)} , (15) 

and 

u{x,z) = x 2 [DiM(r + 2,4,x) + D 2 U(t + 2, 4, x)] exp (-x/2) 

[A sin(ujz) + B cos(luz)} . (16) 

2.2 Bounded plasmas (7 = 0) and C3 = 

Equation Q then reduces to 



x 2 P" + 2nxP' 



n(n — 1) — tx — -x 2 



0. (17) 



4 

This equation can be solved by a procedure similar to that in Sec. 2.1, i.e the 
substitution P = x k W(x), where k is again a root of (|T2|). leads to 

xW» + 2(n + k)W'-(\x + T)w = 0. (18) 

Equation (fTHj) has solutions of the form W(x) = exp(x/2)g(x). For the two roots 
of (fT2"j) . £4 = — n and fc 2 = — (^ + 1) respectively, the solutions are 

W(x) = exp(x/2) [D^i-r, 0, -x) + D 2 U(-t, 0, -a;)] , 

and 

W(x) = exp(a;/2) [D X M{-{1 + r), -2, -x) + L> 2 C/(-(l + r), -2, -x)] . 
The respective solutions for u(x, z) are 

u(x, z) = exp(x/2) [DxM(-r, 0, -x) + D 2 U{-r, 0, -x)] T(z) (19) 
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and 



u(x, z) = exp(x/2)x~ 1 [D x M{-{t + 1), -2, -x) + D 2 U(-(r + 1), -2, -x)} T{z) 

(20) 

with T(z) as given by (fTTH) . 

As in Sec. 2.1, for special values of r and n the solutions of (fTTjl can be 
expressed in terms of simpler classical functions or polynomials, e.g. for n = 
they are expressed in terms of Coulomb wave functions. In this case respective 
static solutions (M p = 0) were constructed in Refs. jHj and [TOj . 

2.3 Bounded plasmas (j = 0) and C3 7^ 

Changing independent variable by 



-2/3 



1 



V = C 3 - ^ 

(jSJ) for n = r = is transformed into the Airy equation: 

d 2 P(v) 

-^U.- v P(ri) = 0. (21) 

The general solution of (|2*T|) is 

P(r ? ) = 4[5 l (r / ) + rf 1 A i (r / )] (22) 

where Ai and i?j are the Airy functions of first and second kind, respectively, and 
d\ and d>2 are constants. An equilibrium symmetric with respect to the mid-plane 
z = then is described by 

u — d 2 [Bi(rj) + d l A i (r])] cos(ciz). (23) 

For C3 = 0, (}2*3*j) takes the simpler form 

u = d 2 



sin(-) +di cos(|) 



cos(Ci2;j. 



In next section solution (|23|) will be employed to evaluate the flow impact on the 
Shafranov shift. 
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3. Flow effects on the Shafranov shift 



In connection with solution (J2HJ) we now specify the plasma container to have 
rectangular cross-section of dimensions a and b and its geometric center to be 
located at R (Fig. 1). Introducing the dimensionless quantities £ = R/Ro, 
( = z/Rq, X = ciR , C = c 2 -Rq? an d 

H = c 3 ~ 2/3 (2c 3 Ce - \ 

(note that c 3 is dimensionless), we require that u vanishes on the plasma bound- 
ary, viz. 

u(H = H ± )=u(( = ±^)=0, 

Ho 

where 

b 



H ± = eg 2/3 (2c 3 Ce ± - \) and £ ± = 1 ± 
This requirement yields the eigenvalues 



R 







\ l = ^(ln + l), 1 = 0,1,2,...; ('' = ('',. A-= 1.2. 



a V 2 

where can be determined by the equations 



A(H+) A(H k ) 

The corresponding eigenfunctions normalized to a reference value u c are given by 

fi w = — = [Bi(H k (£) + D k A t (H k (0)] cos(AiC). 

The simplest eigenfunction corresponding to (k — 1, Z = 0), 

- = mofcC) = [^(gi(0) + £>i^(gi(Q)]cos(AoC) , , 

««(£ = 6i,C = 0) 5 i (i? 1 (£ a ))+ J D 1 A(^l(£a)) ' 1 J 

describes a configuration with a single magnetic axis located on ( a = and 
£ = £ a = 1 + A£ with £ satisfying the equation 

^ +jDl # -°- (25) 

Here, A£ is the Shafranov shift. Contours of constant u and the Shafranov shift 
in connection with (|24|) are illustrated in Fig. 2. 
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For C3 = for which the functions Ai and Bi reduce to cos and sin, respectively, 
the quantities Cf;, and uiq take the simpler forms: = Rokir/^Ab), = 
— tan(Cfc^i) = — tan(Cfc£i) and 



u w = — 

U c 



sm(d£ 2 ) + D 1 cos(Ci£ 2 )J cos(A C)- (26) 
The magnetic axis of the configuration described by is located at (£ = 0, £ = 



l + y/l + P/B®. 

The toroidal current density is given by 



A*ip _ 1 
W ~ / \ 1/2 



R(l- Ml 





|Vw 


2 


dM 2 v 




1 - Ml) 


du 



U=IL 



For Mach numbers of the form M 2 = u m , where m > 1 the profile of j^, is peaked 
on the magnetic axis and vanishes on the boundary. 

On the basis of equation (|25jl . A£ has been determined numerically as a 
function of the flow parameter C3 (Fig. 3). It is recalled that C3 is related to the 
density and the electric field and their variation perpendicular to the magnetic 
surfaces by (0). As can be seen in Fig. 3, A£ is a decreasing function of C3 which 
goes down to zero sharply as c 3 approaches a positive value, this value being 
larger the smaller is the aspect ration Ro/b of the configuration. This result is 
independent of the plasma elongation a/b. It is finally noted that suppression of 
the Shafranov shift by a properly shaped toroidal rotation profile was reported 
inRef. jH|. 

6. Conclusions 

We have constructed exact solutions of the equation describing the MHD 
equilibrium states of an axisymmetric magnetically confined plasma with incom- 
pressible flows [Eq. (HJ)] for S = —ku, where S = d/du [g(d&/du) 2 ], corresponding 
to flows of arbitrary direction. The solutions are based on the form (JJJ) which is 
convenient for applying boundary conditions associated with either unbounded 
plasmas or bounded ones. For k = the solutions are expressed in terms of 
Kummer functions [(Eqs. (fUj) and (|TB|l for unbounded plasmas; (fT^j) and (pU|) 
for bounded ones] while for k 7^ they are expressed in terms of Airy functions 
[Eq. m\- 
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Solution ()23|) has then been employed to study a tokamak configuration of 
arbitrary aspect ratio being contained within a perfectly conducting boundary of 
rectangular cross-section and toroidal current density profile which can be peaked 
on the magnetic axis and vanish on the boundary. In this case it turns out that 
the Shafranov shift is a decreasing function of k which can vanish for a positive 
value k = k c , with k c being larger the smaller the aspect ratio of the configuration 
is. These results demonstrate that the shape of the density and the electric field 
profiles associated with flow and their variation perpendicular to the magnetic 
surfaces can result in a strong variation of the Shafranov shift. 
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Figure captions 



Fig. 1 Boundary with rectangular cross-section in connection with the equilib- 
rium solution (122(1 . 

Fig. 2. The contours illustrate magnetic surface cross-sections in connection with 
solution (|24|) for two different values of the flow parameter C3 which is related to 
the density and electric field and their variations perpendicular to the magnetic 
surfaces. The positions £0 = 1 and £ a refer to the geometric center and the 
magnetic axis, respectively. For C3 = 0.025 the Shafranov shift A£ = (R a — 
Rq)/Rq = £ a — £ m Fig- 2b is drastically reduced in comparison with that for 
c 3 = 0.005 in Fig. 2a. 

Fig. 3. The Shafranov shift A£ = (R a — R )/Ro versus the flow parameter C3 
for various values of the aspect ratio Ro/b. A£ is determined numerically on the 
basis of in connection with the equilibrium solution (|24j) . 
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Figure 1: Boundary with rectangular cross-section in connection with the equi- 
librium solution (|22p. 
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Figure 2: The contours illustrate magnetic surface cross-sections in connection 
with solution (J24j) for two different values of the flow parameter C3 which is related 
to the density and electric field and their variations perpendicular to the magnetic 
surfaces. The positions £ = 1 and £ a refer to the geometric center and the 
magnetic axis, respectively. For C3 = 0.025 the Shafranov shift A£ = (R a — 
Ro)/Ro = £ a — £0 m Fig. 2b is drastically reduced in comparison with that for 
c 3 = 0.005 in Fig. 2a. 
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Figure 3: The Shafranov shift A£ = (R a — Ro)/R Q versus the flow parameter C3 
for various values of the aspect ratio Ro/b. A£ is determined numerically on the 
basis of (J25j) in connection with the equilibrium solution (|24p. 
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